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We present a new approach in considering and including both the perturbative and the 
nonperturbative contributions to the multiplicity ratio r of gluon and quark jets. The new 
method is motivated by recent developments in timelike small-x resummation obtained in 
the MS factorization scheme. A global analysis to fit the available data is also presented. 

1 Introduction 

The gluon-quark multiplicity ratio is defined as r = N g /N q , where N g r q -j is the number of 
hadrons produced in a gluon(quark) jet. A purely perturbative and analytical prediction has 
been achieved by a solution to the equations for the generating functionals in the modified 
leading logarithmic approximation (MLLA) in Ref. PQ up to the so called N 3 LOr in the expansion 
parameter 70 = y/2N c a s /ir i.e. 7$. Here the theoretical prediction is about 10% higher than 
the data at the scale of the Z° vector boson and the difference with the data becomes even 
larger at lower scales. Among the many attempts to predict r numerically, the most successfull 
refers to numerical sulutions to the coupled system of equations of the generating functionals 
for the quark [Zp) and the gluon (Zq) in the MLLA framework (see e.g. [2]). These numerical 
solutions describe well the data only above at relatively high energies [3J 2J [5] . This shows 
that the slope of the multiplicity ratio predicted by this approach tends to be smaller than 
its experimental value. An alternative approach was given in Ref. 6 where equations for the 
derivative of the ratio of the multiplicities are obtained in the MLLA within the framework 
of the colour dipole model. There a constant of integration which encodes nonperturbative 
contributions is fixed by the data. Here a new approach is presented. 

2 The multiplicity ratio in the effective-u; approach 

We consider the standard Mellin-space moments of the coupled gluon-singlet system whose 
evolution in the scale fj, 2 is governed in QCD by the DGLAP equations: 

% ( X ) = ( PZ Pgg ) ( D' B ) ■ (1) 
The timelike splitting functions Py can be computed perturbatively in the strong coupling 
constant: 

(2) 
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where us — N — 1 with N the usual Mellin conjugate variable to the fraction of longitudinal 
momentum x. The functions P^\uj) with k = 0,1,2 appearing in Eq.® in the MS scheme 
can be found in Ref . [5] up to NNLO and in Ref.flO] the NNLL contributions up to 0{a] 6 ) 
in the same scheme. Fully analytical resummed results in a closed form in the MS scheme are 
known at NLL for the eigenvalues of the siglet-gluon matrix jTOJ E] • 

It would be desiderable to fully diagonalize Eq.([T|). However in general this is not possible 
because the contributions to the splitting function matrix do not commute at different orders. 
One is hence enforced to write a series expansion about the LO which in turn can be diago- 
nalized. Therefore we start choosing a basis where the LO is diagonal (see e.g. [12]) with the 
timelike splitting function matrix taking the form: 



P(u;) 



P_ 



where by definition 



3(0), 



P_ 



(3) 



(4) 



and where P±±(uj) are the eigenvalues of the LO splitting matrix. 

Now relating the (g, s) basis to this new (+, — ) basis, we can decompose the singlet and the 
gluon fragmentation function symbolically in the following way: 

£> a (u;,/i 2 ) =D+(w, I i 2 )+D-(w, I jl 2 ); a = s,g. 

According to Eq.([T|) the plus and minus components have the form 



(5) 



u s {n 2 y 



p(0) 



(0) 



where the normalization factors D^{uj,fi 2 ) satisfy 



Mo) 



D g (w, nl) 



D s (w,//o), 



(7) 



with 



(0), 



(0), 



pL°l(u,)-pmu>) 



(0), 



(8) 



The perturbative functions H^(u, /i 2 ) in Eq.© up to NNLO may be represented as 



H±(u:^ 2 ) = 1 



a s (fJ- 2 ) 

47T 

a s (n 2 ) 



4tt 



7(2) 



7(2) 



.(«) 



(9) 



where the functions a , Z. 



(i) 



(2) 

±±,a 



and Z 



(2) 

±=F,a 



with a = q, s in terms of the timelike 



splitting functions up to NNLO in the (+, — ) basis are given by: 



7 (i) 

J ±±,s 



)(0), 



.01 



(10) 



2/3 + Pi°2H-P^M 



(o). 



(i) 

±T,9 



pi°V)-p^H 



(0), 
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Figure 1: 



pert 



J (Q 2 ) with k 



0, 1,2 



(blue, red, green) according to Eq.| |15II , The bands cor- 
respond to rij = 4, 5. 



Figure 2: The multilicity ratio r according to 
Eq.|Q3J with K(n f = 4) = 0.9242 ± 0.0186 (orange) 
and K(n f = 5) = 0.9612 ± 0.0191 (dashed grey). 



and 



m 

J ±±,s 



M = 



J ±T,9 



7 Po 



i=± -I >.=+ 



i=± 



i=± 



i=± 



(ii) 



It is a well known fact that the multiplicity can be obtained from the DGLAP evolution 
equations Eq.([T]) once one is able to take its first Mellin moment ui = N — 1 = 0. This is 
not possible using a fixed order computation because of the presence of singularitues at to = 
due to multiple soft emissions. Resummation of these divergences has been shown to be the 
appropriate thing to do to avoid this problem. This has been shown a long time ago in |13| 
at leading logarithmic accuracy (LL). The algebraic relations in Ref. [H] show that the first 
Mellin moment of the resummed leading logarithmic splitting function Pff (w) can be obtained 
by taking the LO leading singular term and assign an effective value to co : 

P^> = 0) = ^; ^ /= 2P^=0) = J2^ =L382V _ (12) 

Our approach consists in adopting the same procedure also to fix P+^ l (lo = 0) and 
pNNLL^j _ rp ne f ormer quantity is analytically known [TU1 [TT], while the last one is 
known up to the 16th order [10] . In this case we have obtained a numerical estimation of the 
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first Mellin moment of P+^ ll (lu) performing a numerical extrapolation. Our result is: 

Uefj LL = 1-3820 V^T + (0.0059 n/ + 0.8754)a s + (0.0300 n f + 1.0881)a^ 2 , (13) 

which is valid for nt — 4, 5 number of active flavors. 

Neglecting the evolution of the minus component in Eq.© and using Eq.([7]) we arrive at 
our definition of the gluon-quark multiplocity ratio which is given by 

^"<<*-§§^ (14) 

where 

by use of Eq.©. FigQ] shows our results for r^ e L rt (Q 2 ) r^ t L (Q 2 ) and r™ LL (Q 2 ) for n f = 4, 5 
and FigJ5] shows our 90% C.L. fit of K in Eas. (|14ll5j) using the NNLL result for r pert . In our 
analysis we have used the first three terms of the u> expansion for the splitting functions and 
the double counted terms due to resummation have been subtracted. The running of a s has 
been evaluated at NNLO with n/ = 5 and with a s (Mz) = 0.118. The data are taken from the 
summary tables of |15] and references therein and from [16] . 
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